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We reconsider the recently proposed connection between density of states in the so-called "non- 
hcrmitian quantum mechanics" and the localization length for a particle moving in random potential. 
We argue that it is indeed possible to find the localization length from the density of states of a 
non-hermitian random "Hamiltonian" . However, finding the density of states of a non-hermitian 
random "Hamiltonian" remains an open problem, contrary to previous findings in the literature. 
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I. INTRODUCTION 

There are situations in physics in which observables 
can be obtained from properties of non-hermitian oper- 
ators. In this context, the recent work of Hatano and 
Nelson on random "Hamiltonians" with an imaginary 
vector potential has caused considerable interest in so- 
called "non-hermitian quantum mechanics" P-fl^l . Moti- 
vated by the study of the pinning of vortices by columnar 
defects in a superconductor, attention has focused on two 
main questions: What is the spectrum of eigenvalues of a 
non-hermitian "Hamiltonian" , and are the corresponding 
eigenfunctions localized or extended in space? 

In the model introduced by Hatano and Nelson, par- 
ticles are hopping on a lattice with a non-hermitian dy- 
namics governed by the "Hamiltonian" ]jj 

Hh = ICh + ^2 Wj i/jjip^ , (la) 
j 

** = -!E ( eh '^>j+a + ^'^jfa^i ) • ( lb ) 

Here, ipj creates the state at lattice site j, a is a directed 
nearest-neighbor vector, t is the bandwidth, and u>j is 
the random (real) on-site potential. Periodic boundary 
conditions are assumed. The "time evolution" induced 
by "Hh is non-unitary because of the imaginary vector 
potential ih. Numerical simulations by Hatano and Nel- 
son support their conjecture that in the thermodynamic 
limit, the spectrum of the non-hermitian operator ([j]) is 
concentrated on the real axis for energies Ree > £o and 
extends into the complex plane near the center of the 
band, — £o < Ree < £q. They showed that the eigen- 
states with real eigenvalues in the region |Ree| > £o 
are localized, while the eigenstates corresponding to com- 
plex eigenvalues are extended in space |Q. The picture 
that emerges from their analysis is that the energy £q, 
which separates the real and complex eigenvalues, serves 
as a "mobility edge" for the non-hermitian problem (bar- 



ring some unforseen "conspiracy" in which all extended 
eigenfunctions in some energy range So — e < |Ree| < £q 
have real energy eigenvalues). As was shown by Hatano 
and Nelson, the value of the imaginary vector potential 
h where the eigenvalues of the non-hermitian "Hamilto- 
nian" start to pop out into the complex plane is related to 
the localization length £(e) of the problem in the absence 
of the imaginary vector potential, 

\h\=((e a )- 1 . (2) 

It was recently suggested by Hatano ]l0| and Gurarie 
and Zee ]ll[ that the relationship (||) can be inverted, to 
use it as a method to extract the localization length of 
the hcrmitian problem from the support of the spectrum 
of the non-hermitian problem. In this way, knowledge of 
the support of the Density of States (DoS) of the random 
Hamiltonian (Q) as a function of the imaginary vector po- 
tential permits the calculation of the localization length 
as a function of energy in the hermitian case h = |h| = 0. 

To our knowledge, the (ensemble averaged) non- 
hermitian DoS for the random operator (|l|) is known 
only in zero @ and one dimension |5]j7]]. In this paper 
we consider the DoS and its relation to the localization 
length for the so-called Lloyd model Jll|, in which the 
random potentials w$ are independently distributed with 
the Cauchy distribution 

P(w) = I „ 1 „ . (3) 

It is believed that the choice of the Cauchy distribution 
(||) (instead of, say, a Gaussian one) does not modify the 
universal properties of Anderson metal-insulator transi- 
tion in dimension d > 2 [|l4| . 

As shown by Lloyd [O], the advantage of this choice 
of the probability distribution is that the ensemble aver- 
aged DoS of the hermitian problem can be found exactly 
in any dimension d. It has been proposed in the litera- 
ture that the DoS can also be obtained for arbitrary d in 
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the case of the non-hermitian Lloyd model [pj . Accord- 
ing to the relation (0), such a result would permit us to 
find the localization length of the hermitian Lloyd model 
for arbitrary d. We show in this paper that the calcu- 
lation in Q does not give the correct DoS when applied 
to dimensions d > 1. We also discuss the difference be- 
tween the Lloyd model in d — 1 and d > 1 and illustrate 
for the one-dimensional Lloyd model how the localiza- 
tion length for the hermitian problem can be extracted 
from the support of the non-hermitian DoS. Calculation 
of the localization length in the Lloyd model for d > 2, 
however, remains an open problem. 



II. LOCALIZATION LENGTH FROM 
NON-HERMITIAN DOS 

We first discuss how one arrives at the relation (||) be- 
tween the mobility edge eo, the imaginary vector poten- 
tial h, and the localization length £(e) of the hermitian 
problem. 

Hereto we consider, for a given realization of the ran- 
dom potential Wj, a non-degenerate eigenvalue e of the 
hermitian operator Ho with periodic boundary condi- 
tions. (The subscript indicates that the imaginary vec- 
tor potential h is set to zero, i.e., that the Hamiltonian 
Ho is hermitian.) Following Ref. 0], we assume that the 
corresponding eigenstate V?o(j) is localized by the ran- 
dom potential wj, i.e., $o(j) is maximum at a site m and 
decays exponentially far away from m: 



$o(j) ~ exp 



ml 



(4) 



By definition, the exponential decay length £(e) in Eq. 
(f|) is the localization length. Let us now switch on an 
imaginary vector potential ih. As long as h is sufficiently 
small, the wave function 



*(j) - e h 'j $ (j) 



(5) 



is a very good approximation to the exact eigenfunction 
®h(j) °f which adiabatically evolves from <&o(j) as 
h is increased. Although <3> satisfies Hh^ = it is 
not an exact eigenfunction, because it violates the peri- 
odic boundary conditions. The error that one makes is 
of order exp{[|h| — l/£(e)]Z/}. 
Hence, as long as 



|h| < 



1 



(6) 



the wave function 'J will be a good approximation, and 
its energy e will remain real and unshifted (up to an ex- 
ponentially small correction, in principle). 

When the magnitude of the imaginary vector poten- 
tial is larger than the inverse localization length l/£(e), 



the wavefunction (Q) will no longer be a good approxi- 
mation. Both the eigenvalue and eigenfunction undergo 
a qualitative change reflecting the non-hermiticity of the 
"Hamiltonian" . Hence, in the limit of an infinite sys- 
tem size, at |h| = l/£(e), a generic eigenvalue e enters 
the complex plane with unit probability, resulting in the 
relation (|^). 

To justify the inversion of Eq. (Q) to find the localiza- 
tion length £(£•) from the support of the spectrum of Hh, 
we note that for |h| s» l/£(e) eigenfunctions are strongly 
sensitive to the boundary conditions. This sensitivity to 
the boundary conditions causes the phenomenon of level 
attraction M with complex eigenvalues coalescing along 
curves in d = 1, or in compact sets in d > 2 as the 
system size increases. The support of the DoS of Hh ap- 
pears to be self-averaging in the thermodynamic limit, 
i.e., subject to decreasing fluctuations as the system size 
increases. Therefore, the mobility edge eo is well-defined 
for the non-hermitian problem. It has thus been pro- 
posed in jit)) and to relate the real part, eo, of the 
energy at which the first complex eigenvalue appears in 
the spectrum of Hh to the localization length defined by 
Eq. (|). 



III. SINGLE PARTICLE GREEN FUNCTION 
FOR THE NON-HERMITIAN "HAMILTONIAN" 

The advantage of the Cauchy distributed disorder is 
that it allows the exact calculation of the (ensemble av- 
eraged) single particle Green function. In this section 
we discuss whether a similar property exists for a non- 
hermitian system with Cauchy disorder. 

The DoS of the non-hermitian "Hamiltonian" ([!]) is 
computed from the Green function or resolvent 



(7) 



where N is the total number of lattice sites (periodic 
boundary conditions are assumed). The DoS phiz) in 
the complex plane reads M 



7T OZ* 



(8) 



For a hermitian system, when Go (z) is analytic for Im z ^ 
Eq. (||) reproduces the usual DoS concentrated on the 
real axis. 

Let us now consider the ensemble average of the Green 
function Go and the DoS po- Lloyd |D| has shown that 
the average Green function (Go) of the hermitian Hamil- 
tonian Ho is related to the Green function K of the 
non-random Hamiltonian /C [ see Eq. (|l])], 



± 



(9a) 



(9b) 
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The angular brackets denote an average over the random 
disorder potential Wj, 7 is the width of the distribution 
of W} [see Eq. (g)], and 6(x) = 1 (0) for x > (x < 0). It 
follows that the average DoS can be expressed in terms 
of the non-random operator /Co only: 

(poOO) = ^~ [Ko(z - ii) - Ko(z + 17)] SQmz). (10) 

Does the Green function relation Eq. (g) also hold for 
the Green function Gh of the non-hermitian "Hamilto- 
nian" 7i^? The answer is positive, provided 



|Imz| > A, 



(11) 



where A is the imaginary part of the eigenvalue of the 
non-hermitian "Hamiltonian" Tih with the largest imag- 
inary part, 



A = sup {maxfclmefe} 
= maxfclm e' k 

= tsinh(|h||a|) for N > 1. 



(12) 



Here e k (e'k), k = 1, . . . , N are the N eigenvalues of Tlh 
(ICh) for a given realization of the disorder and the supre- 
mum in the first line is taken with respect to all possible 
disorder realizations. 

To see why this is so, we choose to express the Green 
function in terms of a replicated bosonic path integral 



G h (z) = ^Tr J V[<t> c 



1 e J 



(13) 



Here, a is a replica index. The sign in the exponent is 
fixed by the condition that the path integral be conver- 
gent. It is + if Imz > A and — if Imz < — A. If neither 
of these two inequalities holds, i.e., if z lies inside the 
strip Imz| < A, the path integral ( |l3|) cannot be con- 
structed. Averaging Eq. (|l3j ) over the random potential 
Wj is easily done with the Cauchy distribution of Eq. (g) 
if I Imz I > A. In that case, the replicated integrand sat- 
isfies the condition of applicability of Cauchy Theorem 
after closing the contour of integration over w either in 
the upper half plane or lower half plane, depending on 
whether the sign + or — is chosen in Eq. (|l3|). It is thus 
found that, for \lmz\ > A, 

{G h (z)) =J2K h (z±i 1 )9(±lmz), \hnz\ > A, (14a) 
± 

K h (z) = ^—^—. (14b) 
N z - ICh 

As in the hermitian case, the right hand side is expressed 
solely in terms of the non-random resolvent K-h(z). 

Equation ( |l4| ) first appeared in Ref. [pi but without the 
restriction |Imz| > A. The authors of Ref. || applied Eq. 
( |l4| ) to the strip |Imz| < A to obtain the non-hermitian 
DoS in the complex plane 



(ph(z))= o h (z + ij)8(lmz) + a h (z - ij)8(-lmz) 

+ -L [K h (z - i 7 ) -K h (z + i 7 )] 5Qmz), (15) 

where Uh(z) — ^ X d z -Kh(z) is the DoS of the non- 
hermitian problem in the absence of disorder. The DoS 
( |l5| ) corresponds to a non-hermitian DoS coalescing both 
on the real axis (second line) and on a compact set in the 
complex plane (first line). 

The analytical continuation of Eq. to the strip 
|Imz| < A in order to find the DoS can be problematic 
since the Green function Gh is a nonanalytic function of 
z where the DoS is nonzero [compare with Eq. (g)]. It 
can only be justified in the thermodynamic limit in one 
dimension where the non-hermitian spectrum collapses 
to a Id curve. We return to this case in the next sec- 
tion. In all other cases Eq. ( |l5| ) is incorrect. To illustrate 
where it might lead to, we consider Eq. (jlj) in the ther- 
modynamic limit TV — > 00 and extract the length scale 
/(e) = X/h hy locating the edge e(h) at which the wings 
of the spectrum fork into the complex energy plane. Us- 
ing the arguments leading to Eq. (Q), one would identify 
1(e) with the localization length. The dependence of I on 
energy e and dimensionality d is then given by 



cosh 



1(e) 



[A (e) + A 2 (e)] 



A n (e) = y/(\e\-(d- 



n)t) 2 + 7 2 . 



\e\ > (d-l)t, 
\e\ < (d-l)t, 

(16) 



Note that the length scale 1(e) given by Eq. ( |l6| ) is finite 
for all energies and all dimensions. Hence, if Eq. ( p"5| ) 
were true, one would conclude that, irrespective of di- 
mensionality, all states are localized in the Lloyd model 
|nj. This conclusion is not surprising in Id or 2d. In 
fact, the length 1(e) agrees with the localization length 
of the Lloyd model in Id |l5|-[l7|l. In 2d, however, 1(e) 
is much smaller than the weak disorder estimate for a 
Gaussian disorder [P^| . Moreover, in d > 2, such a con- 
clusion contradicts the belief that the existence of large 
tails in the Cauchy distribution does not modify the uni- 
versal properties of Anderson metal-insulator transition 
H|. We return to the issue of dimensions d > 2 and 
the interpretation of 1(e) in section The reason why 
the length scale 1(e) cannot be interpreted as the local- 
ization length for d > 1 is that analytical continuation of 
Eq. ( O) to the strip Im z < A is in general invalid unless 
the DoS is supported on a one dimensional curve, like 
in Id or in the hermitian case h = 0. In particular, we 
conclude that Eq. (|l5|) does not yield the average DoS of 
the non-hermitian extension of Lloyd model in d > 1. 

IV. NON-HERMITIAN DOS FOR ONE CHAIN 

In view of the unreliability of analytic continuation of 
Eq. (|l4j) , it is important to compare Eq. ( |l4| ) with what is 
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known about the spectral properties of the non-hermitian 
"Hamiltonian" TLh from other methods. 

First, we note that in any dimension, analytic contin- 
uation of Eq. (|lj) is certainly wrong in a system of finite 
size: For any finite system and for any dimension the sup- 
port of the averaged DoS (ph(z)) occupies the entire strip 
in the complex energy plane that is excluded in Eq. (11). 
To see this, choose the realization w\ = ■ ■ ■ = wn = V 
with V an arbitrary real number. Equation Jl4|), how- 
ever, results in a DoS with a significantly smaller support: 
it is the DoS of the system without disorder shifted by 
an amount ±7 towards the real axis || . 

What about the DoS in the thermodynamic limit 
N — > 00? Let us first discuss the one-dimensional Lloyd 
model. A discussion of the case d > 1 is postponed to 
the next section. In one dimension, several independent 
approaches have been taken in the literature | 
the Lloyd model, Goldsheid and Khoruzhenko 
shown that the support of the spectrum of Hh is self- 
averaging in the thermodynamic limit and found a DoS 
that coincides with Eq. (|l5|). Hence, the DoS obtained 
from Eq. (|l|) is correct in an infinite one-dimensional 
system, despite the flaws in its derivation. Starting from 
this non-hermitian DoS, one can use the arguments of 
section [H] to identify /(e) with the localization length £(e) 
of the Lloyd model in one dimension. 

We find it instructive to present an alternative deriva- 
tion of Eq. (]l5|) for weak disorder, using the approach of 
Rcf . H , where the support of the DoS was calculated for 
weak Gaussian disorder. In this approach, knowledge of 
the localization length is required to calculate the non- 
hermitian DoS. In the absence of disorder, the energy 
spectrum is parameterized in terms of the (complex val- 
ued) wavenumbers s = 2ir/N J r ih 1 ■ • ■ , 2-K+ih of the plane 
wave states diagonalizing Kh, 




e'(s) — —tcoss. 



(17) 



In Ref. ||, a transfer-matrix approach was used to cal- 
culate the spectrum of Tih for weak non-hermiticity and 
weak disorder to leading order in 1/N. Weak non- 
hermiticity means |h||a| -C 1, whereas weak disorder 
amounts to |sinRes|£ 3> 1, where £ is the localization 
length of the Id hermitian Lloyd model, see Eq. (^6|). To 
leading order in 1/N, it was found that || 



iim S |= \h\ -r 1 

= \h\ 



V^ 2 - (Ree) 2 



(18a) 

0( 7 2 A 2 ). (18b) 



Here, we have expanded Eq. (|16|) to lowest non-trivial 
order in £ _1 and j/t. Hence, we find that the support of 
the spectrum is given by 



in agreement with Eq. ( |l5| ) when d — 1. The variance of 
Ime vanishes like iV^ 1 in the thermodynamic limit. Note 
that Eq. (|l^) relies on very general properties of one- 
dimensional disordered systems through Eq. (18a) and 
only on the specificities of the hermitian Lloyd model 



through Eq. (18t). One dimension is very special in that 
localization length and DoS are closely related jl6| . The 
situation is more intricate in higher dimensions where a 
new length scale, the mean free path, appears besides the 
localization length. 

For small 7, the DoS in the Lloyd model is qualitatively 
different from the non-hermitian DoS in the presence of 
Gaussian distributed disorder Hf7|: In the latter case, 
the uniform shift 7 in Eq. (fL9h should be replaced by the 



non-uniform and much smaller shift 7 2 / ^Jt 2 — (Res) 2 , 
where 7 2 is the variance of the Gaussian distribution. 
To explain the difference, we compare the 7-dependences 
of the localization lengths at the center of the band for 
both disorder distributions. For the Cauchy distribution 
we have £(0) = i/7, while £(0) = t 2 /^ 2 for a Gaussian 
distribution |jl9). Using the relation ( |l^ ) between local- 
ization length and non-hermitian DoS, such a difference 
is directly carried over to the DoS. In fact, the mechanism 
of localization in Id is very different for Cauchy disorder 
compared to that of Gaussian disorder: For a Cauchy 
distribution, the dependence on the disorder strength 7 
of the inverse localization length ~ 7/i follows im- 
mediately from estimating the probability that the dis- 
order potential Wj at an arbitrary site be larger than 
the band width i, in which case the chain is classically 
broken. Thus, localization in the one dimensional hermi- 
tian Lloyd model is not caused by quantum interferences 
effects, but rather by wave functions accomodating to 
large fluctuations in the disorder by vanishing locally. In 
contrast, in the case of Gaussian distributed disorder, 
localization is entirely due to quantum interference. 



V. HIGHER DIMENSIONS: MEAN FREE PATH 



Im e = — t cos(Re s) sinh(Im s) 
= ±h^t 2 - (Ree) 2 T7, 



(19) 



In section III we have shown that analytical continua- 
tion of Eq. (|J into the strip \lmz\ < A yields a length 
scale /(e) that remains finite irrespective of dimension- 
ality. If this analytical continuation were justified, the 
arguments of section || would allow us to identify /(e) as 
the localization length of Lloyd model. Then the Lloyd 
model would not display a metal-insulator transition ir- 
respective of dimensionality. However, as we have seen, 
in general, Eq. (14) cannot be applied inside the strip 
|Imz| < A. This does not mean that the length scale 
/(e) obtained from Eq. (|l4|) is entirely meaningless. In 
this section, we compare the length scale /(e) defined by 
Eq. ([l6]) with other length scales that appeared in pre- 
vious studies of the Lloyd model |^ ^3|, and find that, 
in the limit of weak disorder, /(e) is to be interpreted as 
the mean free path, rather than the localization length. 
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[For strong disorder there is no distinction between mean 
free path and localization length.] The fact that analyt- 
ical continuation of Eq. ([IJ) yields the correct DoS and 
localization length in one dimension is thus simply un- 
derstood as the fact that localization length and mean 
free path coincide in d = 1 . 

To see what is the correct interpretation of the length 
scale 1(e) defined in Eq. Jl6|), we go back to the work of 
Johnston and Kunz [E0|, who considered the quantity 



l' r (z) 



-r 1 Qa.\Gjj+r(z)\) , 



(20) 



where Q(z) = (z — Wo) 1 - They found an expression for 
Eq. (^) valid at least for sufficiently large |Im z\, at fixed 
r and system size. Johnston and Kunz also found that 
if their expression for l' r (z) is analytically continued to 
z->eon the real axis, and the limits of infinite system 
size and separation r are then taken, one obtains a result 
/'(e) which coincides with Eq. (|l6|). Finally, Johnston 
and Kunz hypothesised that l'(s) should coincide with 
the localization length, which one may define by 



1 



]T(7\ = ~ \ J™! r i™ ln \$h3+r 0) I 



(21) 



[We remark that the localization length £(e) is a self- 
averaging quantity, so that the ensemble average in Eq. 
( pl| ) is not necessary] However, as was noticed by Thou- 
less [^lj, MacKinnon |2^] and by Rodrigues, Pastawski, 
and Wcisz ]23| ], this hypothesis is not correct. Moreover, 
Rodrigues, Pastawski, and Weisz showed that the result 
( |l6| ) for /'(e) coincides with the length scale /"(e) defined 

by 



/"(e) 



= — lim 



In 



lim (Gj,j+r(z)) 



(22) 



This length scale is naturally interpreted as the scale on 
which the phase of wave functions is randomnized, rather 
than the length scale on which the amplitude decays ex- 
ponentially. The length /"(e) can also be interpreted as 
the mean free path of a particle for weak scattering. Un- 
fortunately, to find the true localization length £(e), de- 
fined by Eq. (21), remains an open problem in the Lloyd 
model in d > 1. 

A more intuitive picture of what is going on in higher 
dimensions can be obtained along the lines of the last 
paragraph of section IV. The length scale /(e) of Eq. 
( pf ) behaves like 



:\<(d-l)t, 



(23) 



for all dimensions and weak disorder. This implies that 
the calculation scheme of Ref. ||, i.e., application of Eq. 
( |l4| ) to the strip Imz| < A, is predicated on the same 
mechanism as in Id, namely the removal of a site with 



probability 7/t. In contrast to Id, the removal of a site 
does not preclude propagation in d > 2. Instead, the 
removal of sites leads to the usual impurity scattering, 
and we see from Eq. ( p3| ) that /(e) has the semiclassical 
interpretation of an average length for free propagation 
between two impurities, i.e., a semiclassical mean free 
path. 

If we are to interpret Eq. ( |l6| ) for weak disorder as the 
mean free path, we have no reason to exclude that one 
parameter scaling Q applies to Hh in the case h = 0. 
It is widely believed that propagation becomes diffusive 
for a window of length scales beyond the mean free path 
/(e). The upper length scale for diffusive propagation is 
by definition, the localization length. Thus, for the Lloyd 
model in the weak disorder limit we are lead to expect 
the same localization properties as for weak Gaussian dis- 
tributed disorder: (i) localization at all energies if d = 2, 
(ii) existence of a diffusive regime for a window of energies 
centered around e = if d — 3. This expectation is con- 
firmed by the localization properties pi] of a caricature 
of Lloyd model defined by a random hopping tight bind- 
ing Hamiltonian whereby the hopping amplitude takes 
the value t with probability 1 — 7/t and vanishes with 
probability 7/L 



VI. CONCLUSION 

In this paper we have outlined that the localization 
length of a particle moving in a random potential can 
be obtained from the averaged DoS of a non-hermitian 
particle in a random potential. Non-hermitian DoS have 
only been calculated in zero and one dimensions. The 
calculation of the non-hermitian DoS in closed form in 
more than one dimension thus remains an open problem. 
We expect that the non-hermitian trick should provide 
an alternative to numerical calculations of the localiza- 
tion length relying on transfer matrix approaches in strip 
geometries. 
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